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Abstract 

Motivated by the large difference between the direct CP asymmetries Acp{B~ 
-K^K") and Acp{B^ it^K~), we combine the up-to-date experimental information 
on B ^ ttK, ttK* and pK decays to pursue possible solutions with the nonuniversal 
Z' model. Detailed analyses of the relative impacts of different types of couplings are 
presented in four specific cases. Numerically, we find that the new coupling parameters, 
and with a common nontrivial new weak phase ~ —86°, which are relevant to 
the Z' contributions to the electroweak penguin sector ACg and AC7, are crucial to the 
observed "vri^ puzzle" . Furthermore, they are found to be definitely unequal and opposite 
in sign. We also find that Acp{B~ — > p^K~) can put a strong constraint on the new Z' 
couplings, which implies the Z' contributions to the coefficient of QCD penguins operator 
O3 involving the parameter C,^^ required. 
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1 Introduction 



During the past several years, the observed discrepancies between the experimental measure- 
ments and the theoretical predications within the Standard Model (SM) for several observables 
in B ^ ttK decays, the so-called "tt/T puzzle" [1], have attracted much attention. Exten- 
sive investigations both within the SM [2], [3l HI [5], El [7] , as well as with various specific New 
Physics (NP) scenarios [HI El [ID], have been performed. 

Averaging the recent experimental data from BABAR [TT], Belle [12], CLEO [13] and 
CDF [11] , the Heavy Flavor Averaging Group (HFAG) gives the following up-to-date rsults [T5] 



Acp{B- K-Ti^) = 0.050 ± 0.025 , 

^cp(^° ^ K"7r+) = -0.097 ±0.012, (1) 

from which the difference between direct CP violations in the charged and the neutral modes 

AA = Acp{B- K-n^) - Acp{B^ K~n+) = 0.147 ± 0.028 (2) 

is now established at about 5a level. 

Theoretically, it is generally expected that within the SM, these two CP asymmetries 
Acp{B^ TT+K') and Acp{B:[^ n'^K~) should be approximately equal. For example, 
based on the QCD factorization approach (QCDF) [16], the recent theoretical predictions with 
two different schemes for the end-point divergence are 

Scheme I (Scenario S4) [3] , (3) 

Scheme II {nig = 0.5 MeV) [H] . (4) 

Here, the Scheme I is the way to parameterize the end-point divergence appearing in hard- 
spectator and annihilation corrections, by complex parameters X^ h = dy/y = ln(m;,/A)(l + 
Pa,h^^'^^'"), with pa,h < 1 and unrestricted (t)A,H [H]- The Scheme II, as an alternative to the 
first one, is the way to quote the infrared finite gluon propagator to regulate the divergence. It is 
interesting to note that an infrared finite behavior of gluon propagator are not only obtained by 
solving the well-known Schwinger-Dyson equation [TTlfT^fTU] . but also supported by recent Lat- 
tice QCD simulations [20]. However, both of these two schemes suffer the mismatch of given 
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by Eq. Furthermore, within the framework of perturbative QCD approach (pQCD) 
and the soft-colhnear effective theory (SCET) [22], the theoretical predictions read 



pQCD [5] , (5) 




SCET P . (6) 

Obviously, the present theoretical estimations within the SM are not consistent with the es- 
tablished AA. The mismatch might be due to our current limited understanding of the strong 
dynamics involved in hadronic B decays, but equally also to possible NP effects [231 121] • 

In some well-motivated extensions of the SM, additional U{1)' gauge symmetries and asso- 
ciated Z' gauge boson could arise. Searching for the extra Z' boson is an important mission in 
the experimental programs of Tevatron [23] and LHC [2S] • Performing the constraints on the 
new Z' couplings through low-energy physics, on the other hand, is very imporatnt for the di- 
rect searches and understanding its phenomenology. Theoretically, the flavor changing neutral 
current (FCNC) is forbidden at tree level in the SM. One of the simple extensions is the family 
nonuniversal Z' model, which could be naturally derived in certain string constructions [2Z], Eq 
models [28] and so on. It is interesting to note that the nonuniversal Z' couplings could lead 
to FCNC and new CP- violating effect [29] , which possibly provide a solution to the afore men- 
tioned "vri^ puzzle". With some simplifications of the nonuniversal Z' model and neglecting 
the color-suppressed electroweak (EW) penguins and the annihilation amplitudes, Ref. [9] gets 
four possible solutions 

Al ■■ {e"'"', 0l} = {0.0055, 110°} , Bl : {^""^ 0l} = {0.0098, -97°} , with^^^ = ; 
Alr : {e"'"' = e'''', <Pl} = {0.0104, -70°} , Blr : {^"^^ = 0l} = {0.0186, 83°} . (7) 

However, the corresponding prediction Acp{B~ n^K~) = — 0.03±0.01 of solution A^ and 
Alr in Eq. (I7j) is obviously inconsistent with the up-to-date experimental data 0.050 ± 0.025. 
Moreover, the annihilation amplitudes, which could generate some strong-interaction phases, 
are important for predicting CP violations. 

Based on the above observations, in this paper we shall adopt the QCDF approach and 
reevaluate the effects of the nonuniversal Z' model on these decay modes with the updated 



experimental data. Furthermore, since the B ttK* and pK decays also involve the same 
quark level b — ^ sqq {q = u, d) transition, it is necessary to take into account these decay modes. 

In Section 2, we provide a quick survey of -B ^ -kK, ttK* and pK decays in the SM within 
the QCDF formalism; our numerical results, with two different schemes for the end-point 
divergence, are also presented. In Section 3, after reviewing the nonuniversal Z' model briefly, 
we present our analyses and numerical results in detail. Section 4 contains our conclusions. 
Appendix A recapitulates the decay amplitudes for the twelve decay modes within the SM [3]. 
Appendix B contains the formulas for hard-spectator and annihilation amplitudes with the 
infrared finite gluon propagator [8]. All the theoretical input parameters are summarized in 
Appendix C. 



2 The SM results with two schemes for the end-point 
divergence. 

In the SM, the effective Hamiltonian responsible for b ^ s transitions is given as 



r r ^° 



V2 



i=3 



+h.c., (8) 

where VqbV*g {q = u, c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements [30] , Ci the Wilson coefficients, and Oi the relevant four-quark operators whose explicit 
forms could be found, for example, in Refs. [HTj. 

In recent years, the QCDF approach has been employed extensively to study the hadronic 
B-meson decays. The B — * ttK, tiK* and pK decays have been studied comprehensively within 
the SM in Refs. [2|, [3l HI [32] , and the relevant decay amplitudes within this formalism are shown 
in Appendix A. It is also noted that the framework contains estimates of the hard-spectator and 
annihilation corrections. Even though they are power-suppressed, their strength and associated 
strong-interaction phase are numerically important to evaluate the branching ratio and the CP 
asymmetry. However, unfortunately, the end-point singularities appear in twist-3 spectator and 
annihilation amplitudes. So, how to regulate the end-point divergence becomes important and 
necessary within this formalism. Here we shall adopt the following two schemes: 



Scheme I: Parametrization 

As the most popular way, the end-point divergent integrals are treated as signs of infrared 
sensitive contributions and phenomenologically parameterized by [21 E] 

r*^X. = (l + p.e'*^)lnfi. r.y'^^-hx.f, (9) 

Jo y Afe Jo 2/ 2 

with Ah = 0.5 GeV, Pa ^ 1 and 0^ unrestricted. is treated in the same manner. The 
different choices of and (pA correspond to different scenarios as discussed in Ref. [3], and S4 
is mentioned as the most favorable one. It presents the moderate value of nonuniversal anni- 
hilation phase (pA = —55° (PP), —20° (PV) and —70° (VP). Conservatively, in our calculations 
we quote ±5° as their theoretical uncertainties. Taking pA = ^ and Xa,h universal for all decay 
processes belonging to the same modes (PP, PV or VP), we present our numerical results of 
branching ratios and direct CP asymmetries for B ttK, ttK* and pK decays in the third 
column of Tables [2] and [3l respectively. 

As is known, the mixing-induced CP asymmetry A^p is well suited for testing the SM 
and searching for new physics effects. For example, the investigation of mixing-induced CP 
asymmetries in penguin dominated n^Kg and 5° P^Kg decay modes has attracted 

much attention recently [331 EH ESI ESj- After neglecting the Kq — Kq mixing effect, the 
mixing-induced asymmetry could be written as 

A-(5° ^ /) = , (/ = vr°A^°,p°A^°), (10) 

with Xf = —exp{i aTg[r^^iy^]}Af / Af in our phase convention. Our numerical predictions are 
listed in Table HI which agree with the measurements within large experimental errors. 

Scheme II: Infrared finite dynamical gluon propagator 

In our previous paper [S], we have thoroughly studied the end-point divergence with an in- 
frared finite dynamical gluon propagator. It is interesting to note that recent theoretical and 
phenomenological studies are now accumulating supports for a softer infrared behavior of the 
gluon propagator [THl EZl EH]- Furthermore, the infrared finite dynamical gluon propagator, 
which is shown to be not divergent as fast as ^, has been successfully applied to the hadronic 
B-meson decays [391 HQ] . In our evaluations, we shall quote the gluon propagator derived by 
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Cornwall (in Minkowski space) [T7] 

— M^\q^) + ze 

where q is the gluon momentum. The corresponding strong coupling constant reads 

where /3o = 11 — |n/ is the first coefficient of the beta function, with n/ being the number of 
active quark flavors. The dynamical gluon mass square M^(g^) is obtained as [T7] 



In 



~f? 



In'^ 



12 
■ 11 



(13) 



where is the effective gluon mass and Aqcd = 225 MeV. In Ref. [8], we present our sugges- 
tion, nig = 0.50 ±0.05 GeV, which is a reasonable choice so that most of the observables (except 
for Acp{B 71^ K~)) are in good agreement with the experimental data. In this way, we find 
that the hard-spectator scattering contributions are real, and the annihilation contributions 
are complex with a large imaginary part [8j. Our numerical predictions for branching ratios, 
direct CP asymmetries and mixing-induced CP asymmetries are listed in the fourth column of 
Tables [21 El and HI respectively. 

Although numerically these two schemes have some differences, both of their predictions 
are consistent with most of the experimental data within errors. However, as expected in the 
SM, we again find that Acp(5,7 tx^K^) = -0.041 ± 0.008 (-0.100 ± 0.008), are very close 
to Acp{B^ TT+K-) = -0.077 ± 0.009 (-0.116 ± 0.008) in the first (second) scheme. So, it 
is still hard to accommodate the measured large difference AA in the SM within the QCDF 
formalism, irrespective of adopting which scheme. In the following, we pursue possible solutions 
to this problem with a family nonuniversal Z' model 



3 Solution to the "tt/T puzzle" with nonuniversal Z' model. 
3.1 Formalism of the family nonuniversal Z' model 

A possible heavy Z' boson is predicted in many extensions of the SM, such as grand unified 
theories, superstring theories, and theories with large extra dimensions. The simplest way to 
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extend the SM gauge structure is to include a new U{1) gauge group. A family nonuniversal Z' 
model can lead to FCNC processes even at tree level due to the non-diagonal chiral coupling 
matrix. The formalism of the model has been detailed in Ref. [29]. The relevant studies in the 
context of B physics have also been extensively performed in Refs. [HI SSI SS] . 

After neglecting the Z — Z' mixing with small mixing angle 6 ~ (9(10~^) [S], and taking 
all the fields being the physical eigenstates, the Z' part of the neutral-current Lagrangian can 
be written as [29] 

C = -g'J'^Z'^ , (14) 

where g' is the gauge coupling constant of extra U'{1) group at the EW Miy scale. The Z' 
chiral current is 

j; = ^PaAiBihPL + (5f , (15) 

where is the mass eigenstate of chiral fields and Pl,r = (1 =F 75)/2. The effective chiral Z' 
coupling matrices are given as 

Bf = V,xe,xV^j,, {q = u,d;X = L,R). (16) 

With the assumption of flavor-diagonal right-handed couplings , the Z' part of the effective 
Hamiltonian for b sqq [q = u, d) transitions can be written as [9] 

-H^f = ^(^)'i^f.(^-&)v'-AE + K(^9)v^a) + h.c. , (17) 

where gi = e/ {sin 6w cos^vy) and Mz' the new gauge boson mass. It is noted that the forms of 
the above operators already exist in the SM. As a result, Eq. (fT7|) can be modified as 

n%f = -^VtlVtb 1](AC303^ + AC.OI + + AC,0I) + h.c. , (18) 

where 0^{i = 3, 5, 7, 9) are the effective operators in the SM, and ACi the modifications to the 
corresponding SM Wilson coefficients caused by Z' boson, which are expressed as 



AC, 



9,7 



SVtlVtb 'giMz 
in terms of the model parameters at the Mw scale. 
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^i^rB^AB'^f-Bif), (19) 



Generally, the diagonal elements of the effective coupling matrices B^^^ are real as a result of 



the hermiticity of the effective Hamiltonian. However, the off-diagonal ones of i?^f, can contain 
a new weak phase 0^- Then, conveniently we can represent ACj aJl] 

AC, J = 4^|J^e^^'^^e^^N (20) 
where the real NP parameters (^^^'^^^ ^ll,lr defined, respectively, as 

^LL,LR 



_1 , g'Mz .2 I I , ^n-RL,R. 

<Pl = Arg[i?f,]. (21) 

It is noted that the other SM Wilson coefficients may also receive contributions from the Z' 
boson through renormalization group (RG) evolution. With our assumption that no significant 
RG running effect between and Mw scales, the RG evolution of the modified Wilson 
coefficients is exactly the same as the ones in the SM [31], HI] . For simplicity, we define 

X = ^LL^i4>L ^ Y = S,^^e'^^ . (22) 
The numerical results of Wilson coefficients in the naive dimensional regularization (NDR) 



scheme at the scale /i = mf, (/i/j = ^/KfJr^) are listed in Table [TJ The values at the scale Hh-, 
with TTifo = 4.79 GeV and = 500 MeV, should be used in the calculation of hard-spectator 
and weak annihilation contributions. 



3.2 Numerical analyses and discussions 

With the theoretical formulas and the input parameters summarized in Appendix A, B and 
C, we now present our numerical analyses and discussions. Our analyses are divided into the 
following four cases with different simplifications for our attention, namely, 

• Case I: With the simplifications B^f ~ -25^/ («.e., (■^'^--f'^ = 0) and = 0, 
^ For comparison, we take the same phase convention as Ref. j9j. 
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Table 1: The Wilson coefficients Ci within the SM and with the contribution from Z' boson 
included in NDR scheme at the scale /i = mb and = V-^hmb- 



Wilson 






coefficients 


(jSM 


ACf 




CSM 


ACf 


Ci 


1.075 


-0.006X 




1.166 


-0.008X 


C2 


-0.170 


-0.009X 




-0.336 


-0.014X 


Cs 


0.013 


0.05X - O.Oiy - 2.20X' 


- o.o5y' 


0.025 


O.llX - 0.02y - 2.37X' - 0.12y' 


Ci 


-0.033 


-0.13X + O.Oiy + 0.55X' 


+ o.o2y' 


-0.057 


-0.24X + 0.02y + 0.92X' + 0.09y' 




0.008 


0.03X + O.Oiy - 0.06X' 


- i.83y' 


0.011 


omx + o.o2y - o.iox' + o.o9y' 


Ce 


-0.038 


-o.i5x + o.oiy + 0.1X' 


- o.ey 


-0.076 


-0.32X + 0.04y + 0.16X' - 1.26y' 


C7/ Qem 


-0.015 


4.18X-473y + 0.25X'- 


f i.27y' 


-0.034 


5.7X - 459y + OAX' + 1.7y' 




0.045 


1.18X - 166y + 0.01X' - 


f 0.56y' 


0.089 


3.2X - 355y + 0.2X' + 1.5y' 




-1.119 


-561X + 4.52y-0.8X' 


+ o.4y' 


-1.228 


-611X + 6.7y - 1.2X' + o.6y' 


Clo/aem 


0.190 


118X -0.5y + 0.2X' - 


o.o5y' 


0.356 


207X - 1.4y + 0.5X' - O.IY' 


C77 


-0.297 






0.360 




Csg 


-0.143 






-0.168 





• Case II: With the simplifications fi^f - -25^/ only {i.e., (^^^^^ = 0), 

• Case III: Taking ~ -2i?J^ (i.e., (^^ ~ 0), and leaving (^^ and ^^^'^^ arbitrary, 

• Case IV: Without any simplifications for B^:^ and Bj^f^, i.e., arbitrary values for C,^^^^^ 
and ^^^'^^ are allowed. 

Our fitting is performed with the experimental data varying randomly within their 2(7 error- 
bars, while the theoretical uncertainties are obtained by varying the input parameters within 
the regions specified in Appendix C. In addition, we quote the Scheme II (taking irtg = O.SGeV) 
to regulate the appearing end-point divergences. 

With the assumption B^^^ ^ —2B^^ and neglecting the color-suppressed EW penguins and 
the annihilation amplitudes, four possible solutions Eq. ([71) to the "vrii' puzzle" are obtained 
in Ref. [9]. It is still worth to recheck these solutions with the updated experiment data and 



9 



Table 2: The CP-averaged branching ratios (in units of 10~^) of -B ^ ttK, ttK* and pK decays 
in the SM with two end-point divergence regulation schemes, and in the nonuniversal Z' model 
with four different cases. 

Decay Mode Exp. SM Z' model 

data Scheme I Scheme II Case I Case II Case III Case IV 

B~^7T~K^ 23.1 ±1.0 19.0 ±2.5 23.4 ± 3.9 23.3 ± 0.7 23.3 ± 0.6 23.2 ± 0.6 23.3 ± 0.7 

B-^TT^K" 12.9 ±0.6 10.5 ±1.3 12.7 ±2.0 12.5 ± 0.6 12.6 ± 0.6 12.5 ± 0.5 12.6 ± 0.6 

b'^^^tt+K- 19.4 ±0.6 16.2 ± 2.2 20.1 ±3.4 19.9 ± 0.5 19.8 ± 0.5 19.9 ± 0.5 20.0 ± 0.5 

^ 7r°K° 9.8 ±0.6 7.3 ±1.1 9.3 ± 1.7 9.4 ± 0.6 9.5 ± 0.6 9.1 ±0.4 9.1 ±0.4 

B- ^ TT~K*° 10.0 ±0.8 11.7 ±1.2 10.3 ±3.3 8.4 ±1.0 8.5 ± 0.9 8.7 ±0.6 8.6 ±0.7 

B-^TT^K*- 6.9 ±2.3 7.0 ±0.7 6.0 ± 1.8 4.7 ±0.6 4.7 ±0.5 4.9 ± 0.3 4.8 ± 0.3 

B°a^7T+K*- 10.3 ±1.1 9.9 ±1.1 9.2 ±2.8 7.5 ± 1.0 7.7 ±0.9 8.0 ± 0.6 8.0 ± 0.6 

B°a^7T°K*° 2.4 ±0.7 4.1 ±0.5 3.9 ± 1.3 3.6 ± 0.5 3.7 ±0.4 3.5 ± 0.4 3.5 ± 0.4 

B- ^ p-K° 8.0t\i 5.2 ±0.9 10.6 ± 3.7 9.6 ± 1.4 9.7 ±1.3 10.6 ±1.3 10.7 ±1.5 

B- ^ p^K- 2,MtQ% 2.5 ±0.4 5.4 ± 1.6 4.22 ± 0.62 4.47 ± 0.63 4.7 ±0.6 4.8 ±0.7 

5° ^ p+K- ^.Qtil 6.3 ±1.0 13.0 ±3.8 10.8 ±1.4 10.9 ± 1.4 11.9 ±1.4 12.5 ± 1.8 

B^^P^k'^ bAtil 3.7 ±0.5 7.3 ±2.1 7.1 ±0.9 7.4 ± 0.9 6.8 ± 0.9 6.9 ±1.0 



taken into account the neglected corrections. Furthermore, the possible solutions may also 
suffer strong constraints from B ttK* and pK decays, since they are also mediated by the 
same quark level b —>■ sqq transitions. 

Case I: With the simplifications 5^;^ ~ -2B^f{i.e., (^^'^^ = 0) and = 

In this case, assuming B^^^ ~ —2B^f^ as in Ref. [9], the NP effect primarily manifests itself 
in the EW penguin sector and the Z' contribution to the Wilson coefficients Eq. (1191) can be 
simplified as 

AC3,5 = 0, 

AC,. = 4^e^^'^V^^ withe-'-^=(f^)^|-^|5,Y. (23) 
vts^tb gi^y^z' vtsnb 

As shown in Fig. □ (a), taking = 0.004 and = 0, we find that Acp{B- tt^K') is 
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Table 3: The direct CP asymmetries (in unit of 10 ^) oi B ^ ttK, ttK* and pK decays. The 
other captions are the same as Table. O 

Decay Mode Exp. SM Z' model 



data Scheme I Scheme II Case I Case II Case III Case IV 





TT K 


0.9 ±2.5 


0.4 ±0.1 


0.04 ±0.07 


-1.6 ±0.3 


-2.7 ±0.9 


5.2 ±0.5 


5.1 ±0.6 


Bu- 




5.0 ±2.5 


-4.1 ±0.8 


-10.0 ±0.8 


2.4 ± 1.6 


2.3 ± 1.5 


0.9 ±0.7 


1.2 ±0.9 


Bd - 


TT+K~ 


-9.8t\i 


-7.7 ±0.9 


-11.6 ±0.3 


-11.7±0.3 


-11.0 ±0.7 


-10.5 ± 1.1 


-10.5 ± 1.2 


Bd - 




-1 ± 10 


-1.5 ±0.3 


0.7 ±0.3 


-17± 2 


-18 ±2 


-6± 2 


-6± 2 


Bu - 


-17*^ 


9+6.7 
^"6.1 


0.6 ±0.1 


0.09 ±0.15 


-2.1 ±0.4 


-3.3 ±0.5 


-0.6 ± 2.4 


-3.0 ±6.7 


Bu - 




4 ±29 


-6 ±2 


-37 ±9 


6.8 ± 7.1 


9.1 ±7.2 


-17±4 


-18 ±6 


Bd - 
R° 

Bd - 


-> 7r+K*~ 
-> tt 'K 


-25± 11 
-15± 12 


-13 ±2 
-4± 1 


-43 ± 10 
4±2 


-48 ±3 
-58 ±9 


-46 ±3 
-62 ±9 


-49 ±3 
-34 ±7 


-50 ±5 
-36± 11 


Bu - 


p K 


-12± 17 


0.4 ±0.1 


0.5 ±0.2 


1.5 ±0.1 


-0.15 ±0.7 


5.3 ± 1.1 


6.5 ±4.5 


Bu - 
R° 

Bd - 


p'^K- 
p+K~ 


15 ±6 


57.3 ±5.8 
36 ±4 


42.3 ±9.5 
29 ±6 


-36 ± 10 
31 ±3 


-46 ± 12 
33 ±3 


27 ±4 
25 ±2 


27 ±5 
25 ±2 


R° 

Bd - 


.p^K' 


1 ± 20 


-2.1 ± 1.3 


-2.4 ± 1.4 


45 ±5 


50 ± 5 


8±3 


9±4 



Table 4: The mixing-induced CP asymmetries (in unit of 10 ^) of B^ —>■ tt^K^ and p^Kg 
decays. The other captions are the same as Table. O 



Decay Mode 


Exp. 

data 


SM 

Scheme I Scheme II 


Case I 


Z' 

Case II 


model 
Case III 


Case IV 


5^ - p'K^s 


57± 17 
63tll 


77 ± 2 77 ± 2 
60 ± 2 66 ± 2 


46 ±6 

87 ±2 


44 ±6 
84 ±3 


61 ±3 
85 ±3 


62 ±5 
86 ±9 



enhanced to be consistent with the experimental data when (pL ~ —90°. Moreover, Acp{B^ — >• 
ir'K^) and Acp{B^ tx^K~), which agree roughly with the experimental data in the SM, are 
not sensitive to the parameter . So, a possible solution to the observed ^^ttK puzzle" Eq. ([2D 
in Case I is naively favored. 
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Figure 1: The dependence of Acp{B — > vrK) on the new weak phase (^l for the values of 
^LR^ ^LL^ g^j^i^ ^LR marked by the legends. 



Taking BiB nK) and Acp{B ttK) as constraints on and (pi, the allowed region 
for these two parameters are shown in Fig. [3] and the corresponding numerical results are listed 
in Table. El z.e., = (3.96 ± 0.70) x lO'^ and 0l = -88° ± 7°. Our result confirms that the 
solution Bl in Eq. is helpful to resolve the "vrii' puzzle" (note that a bit of difference might 
be due to the fact that the annihilation corrections are not included in Ref. [9]). However, the 
solution Al is excluded by the updated experimental data Acp{B~ ti^K^) = 0.050 ± 0.025 
as indicated in Fig. [T] (a). 

With = (3.96 ± 0.70) x 10^^ and (pi = —88° ± 7° as input parameters, we present our 
predictions for B{B nK*,pK), Acp{B nK\pK) and A^)?(fi° n^Ks.p^Ks) in the 
fifth column of Tables. [21 Eland HI respectively. We can see that most of them are consistent with 
the experimental data within 2a. Especially, the predicted A^^p{B^ tx^Ks) = 0.46 ± 0.06 
is very close to the measurement 0.57 ± 0.17 [15]. However, the prediction for Acp{B^ — >■ 
p^K^) = —0.36 ± 0.10 presents a large discrepancy (larger than 6a errors) with the current 
experiment data 0.419^0 [IS], which is also shown in Fig. [21 (a). This fact implies that 
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Figure 2: The dependence of Acp{B — > pK) on the new weak phase (j)L- 



-50 



(U 

^-100 



-150 



^L^[xlO"^] 



Case I 



Figure 3: The ahowed regions for the parameters and <?5>l in Case I. 



Acp{B p^K ) can provide a strong constraint on the Z' couplings, at lease in Case I, and 
some more general Z' models might be required to explain all of these measurements. 
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Figure 4: The allowed regions for the parameters ^^^'^^ and (pL in Case II. 



Case II: With the simphfication B^f ~ -25^/ only (z. e. = 0). 

It is interesting to note that, as shown in Fig. [11(b), a region of minus with (j)L 90° can 

bridge the discrepancy of Acp{B^ tt^K~) between theoretical predictions and experimental 
data. Moreover, it is also possible to moderate the problem of Acp{B~ P^K ) induced by 
as shown in Fig. [2] (b). So, in Case II we give up the simplification = and pursue 
possible solutions to these discrepancies. 

Taking B{B — > irK) and Acp{B irK) as constraints, we present the allowed regions 
for and (pi in Fig. HI Unfortunately, we find that the required region of minus C,^^ 

with (pL ~ —90° is excluded by Acp{B^ 7i~^K~), because it will induce a large negative 
Acp{B^ — i> TT^K') as shown in Fig. [T] (b). In addition, as shown in Fig. [T] (b), the region of 
plus C,^^ with (pL ~ —90° is helpless to resolve the "vr/^ puzzle" . The Z' effects are therefore 
still dominated by large and the problem of Acp{B~ P^K~) induced by still exist. 

In fact, with and having the same sign, the corresponding Z' contributions coun- 
teract with each other in the B tx^K~ decay as shown in Figs. [1] (a) and (b). It is also easily 
understood from the expression for the effective coefficient ^y^{P P) = a^{PP) — a^{PP) [5], 
which involves the leading-order Z' contribution in this case. Thus, we conclude that any at- 
tempt to explain the B — > nK anomaly in the non-universal Z' model with the assumption 
= = ^, as made in Ref. [IS], is frangible and excluded in our case. 

In a word, although the Z' contributions with a positive C,^^ or a negative and 0l ~ —90° 
are helpful to bridge the discrepancy of Acp{B^ n^K'), they would induce the unmatched 
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Table 5: The numerical results for the parameters ^^^'^^^ ^ll,lr ^^^^ different 

cases. The dashes mean that the corresponding parameters are neglected in each case. 



Parameters Case I Case II Case III Case IV 





xlO" 




3.96 ±0.70 4.32 ±0.75 


1.52 ±0.24 


1.65 ±0.35 




;xiO' 




— 0.21 ±0.15 


-0.53 ±0.13 


-0.54 ±0.15 


en 


;xiO' 






-11.8 ±3.1 


-14.6 ± 7.1 




;xio- 








1.04 ±2.70 








-88° ± 7° -88° ± 7° 


-86° ± 14° 


-85° ± 16° 



Acp{B- p^K-) and Acp(5° n+R-), respectively. Thus, with both B{B txK) and 
Acp{B irK, pK) as constraints, our results indicate that all of the parameter spaces in Case 
I and Case II are excluded with the assumption -B^l^ — —2B^^. As an alternative, in the 
following, we proceed to pursue possible solutions to these observations by considering the Z' 
contributions to the QCD penguins AC3,5. 

Case III: Taking B^^ ~ -2B^^{i.e., ~ 0), and leaving and ^^^'^^ arbitrary. 

As shown in Fig. [T] (c), we find that the variation trends of Acp{B^ n^K") and Acp{B~ — >■ 
7r^K~) are always the same, indicating that the Z' contributions in this case could not give a 
solution to the observed "vri^' puzzle" directly, as well as the unmatched AcpiB'^ P^K~) in- 
duced by . However, it is interesting to note that, with ~ —90°, both Acp{B^ — > ti^K~) 
and Acp{B^ 7i^K~) could be enhanced simultaneously, which may relax the constraints on 
As mentioned in Case II, a negative is favored by the "vrii' puzzle" and can mod- 
erate the problem of Acp{B^ — > p^K~) induced by So, the parameter C,^^ may play an 
important role. 

With E{B — > ttK), Acp{B ttK) and Acp{B pK) as constraints, the allowed regions 
for C,^^, Q^^ and (pL are shown in Figs. [5l We find that none of S,^^ , and could 
be neglected. Especially, the part moderates the contradictions caused by C,^^ and 
Furthermore, it is interesting to note that our predictions for B{B ttK*,pK), Acp{B — >• 
ttK*) and n^Ks.p^Ks), listed in Tables El [3] and H respectively, are all consistent 
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Figure 5: The allowed regions for the parameters ^^^'^^^ , and 0/, in Case III. 



with the experimental data within 2a. 

Case IV: Without any simplification of B^^^ and i.e., arbitrary values of ^^^'^^ 

and ^^^'^^ are allowed. 

More generally, we give up any assumptions of the couplings 5^^^ and B^. Then, there 
are five arbitrary NP parameters. As in Case III, we take B{B ttK), Acp{B -kK) and 
Acp{B pK) as constraints and present the predictions for the other observables. 

The allowed regions for ^^^-^^^ ^ll,lr shown in Fig. [6|, while the numerical 

results are listed in the last column of Table We find that, similar to Case III, the values of 
^LL,LR g^j,g definitely nonzero. The values of C,^^ is a little larger than the one in Case III, due 
to the interference effect caused by the parameter C,^^. Our predictions for B{B 'KK*,pK), 
Acp{B ttK*) and ^ n^Ks, p'^Ks), listed in Tables [21 E] and H respectively, are 
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Figure 6: The allowed regions for the parameters ^^^-^^^ qLL,lr Case IV. 



consistent with the experimental data within 2a. 

4 Conclusions 

Motivated by the recent observed large difference A A between Acp{B^ ir^K^) and Acp{B^ - 
K^-K^), we have investigated the effect of family non-universal Z' model and pursued possi- 
ble solutions to the observed "vri^' puzzle". Moreover, we have also taken into account the 
constraints from the B — > irK*, pK decays, which also involve the same quark level b — sqq 
{q = u, d) transitions. Our main conclusions are summarized as: 

• The Z' contributions to the coefficients of operators O-j and Og {^^^ and ^^^) with 
0L ~ —86° are crucial to bridge the discrepancy of Acp{B~ tt^K~) between theoretical 
prediction and experimental data. However, they are definitely unequal and opposite in 
sign. 
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• The Z' contributions to the coefficients of QCD penguins operator O3 related to C^i^i^ are 
required to moderate the contradiction of AcpiB' P^K~) and Acp{B^ 7r^K~) 
to thier experimental values induced by and respectively, even though they are 
helpless to resolve the observed "vri^ puzzle" . On the other hand, the Z' contributions to 
C^iC^^) are inessential. 

• For all of the four cases, a new weak phase associated with the chiral Z' couplings, with 
a value about —86°, is always required for the "vrii' puzzle". 

Combing the up-to-date experimental measurements of B ^ ttK, piK* and pK decays, the 
family non- universal Z' model is found to be helpful to resolve the observed "vri^" puzzle". It 
is also reminded that more refined measurements of the mix-induced CP asymmetries in the 
n^Ks and p^Ks decays are required to confirm or refute the NP signals. In the following 
years, the precision of measurements for these observables is expected to be much improved, 
which will then shrink and reveal the Z' parameter spaces. 

Note added: When the paper is finished, we are aware of the interesting paper by Barger 
et al. /5^. Although our topics are very similar, we have taken into account of not only the 
CP asymmetries but also the branching ratios of the correlated decay modes to constrain Z' 
couplings. Moreover, our approaches for the hadronic dynamics are different. 
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Appendix A: decay amplitudes in the SM with QCDF 

The decay amplitudes for B txK decays are recapitulated from Ref. [3] 
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B~-*ttOK- 



p=u,c 



3,EW 



ylSM 



EW 



(25) 
(26) 



p=u,c 



p=u,c 



5p« "2 + 2'^3,EW 



(27) 



where the exphcit expressions for the coefficients af = af (M1M2) and /5f = /5f (M1M2) can also 
be found in Ref. [3]. Note that expressions of the hard-spectator terms Hi appearing in and 
the weak annihilation ones appearing in Pf should be replaced by our recalculated ones listed 
in Appendix B. The decay amplitudes of 5 ^ ttK* and B pK decays could be obtained 
from the above results by replacing (irK) {irK*) and (ttK) — > (pK), respectively. 



Appendix B: The hard-spectator and annihilation correc- 
tions with the infrared finite gluon propagator 

With the infrared finite gluon propagator to cure the end-point divergences, the hard-spectator 
corrections in i? ^ PP and PV decays can be expressed as [H] 



dxdyd^ — - — $Bi(0'^A/2(a;) 
Jo s 



_l_ ^ 



(28) 



for the insertion of operators Qi=i-A; 



9,10, 



Jo ^ 



Ml 



0mi(Z/) 



x{y + uj^{q^)/0 x{y + uj^{q^)/0 

(29) 



for Qi=5j, and Hi(MiM2) = for Qi=6,8- When both Mi and M2 are pseudoscalars, the final 
building blocks for annihilation contributions can be expressed as [Sj 

X 1 



A\ = IT dxdyas{q' 
Jo 

+^ . .. rf rf 0^,(|/)0^,(x) \ , 



+ 



{xy — uj'^{q'^) + ie)(l — xy) {xy — uo'^{q^) + ie)x 



xy — uj'^{q'^) + ie ^ ^ 



(30) 
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{xy — uj'^{q^) + ie){\ — xy) {xy — uj'^{q^) + ie)y 



Ml 



2(n2\ I ,V X X 



xy — uj^{q^) + le 



rfvf 0^,(2/)0^,(x)L (31) 



-rf (x)0„,(x)$M.(y)L (32) 



(xy — Ci;2(g2) + «e)(l — x|/) ^ 

4 = 4 = 0, (33) 



^ ^ 2(1 + x) 



aI = TT dxdyas{q^){j^ "^^7"^, . , rf i0mi(l/)$M,(x) 



+ ^^.Ltf , (x)0^,(x)$M,(y) 1^ . (34) 



/o [ (xy — 07^(^2) + ie)x 

2(1 + ?/) 
[xy — uj'^{q'^) + ie)y ^ 

When Ml is a vector meson and M2 a pseudoscalar, the sign of the second term in A\, the 
first term in Ag, and the second terms in A\ and /I3 need to be changed. When M2 is a vector 
meson and Mi a pseudoscalar, one only has to change the overall sign of A\. 

Appendix C: Theoretical input parameters 
CI. CKM matrix elements 

For the CKM matrix elements, we adopt the Wolfenstein parameterization [l6] and choose the 
four parameters A, A, p and 77 as [17] 

A = 0.798l°:°f7, A = 0.225211™, 'p = Q .Ult^'.^il, = 0.340 ± 0.016, (35) 

with p = p (1 — ^) and f] = rj{l — 



C2. Quark masses and lifetimes 

As for the quark masses, there are two different classes appearing in our calculation. One type 
is the current quark mass which appears in the factor r^^ through the equation of motion for 
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quarks. This type of quark masses is scale dependent and denoted by rriq. Here we take 

frisifi) /Wigifi) = 27.4±0.4 |1H], m,(2 GeV) = 87±6 MeV [IE], mb(mb) = 4.20+[J:J^ GeV 

(36) 

where rnq{ij) = {friu + ffidjifi) /2, and the difference between u and d quark is not distinguished. 

The other one is the pole quark mass appearing in the evaluation of penguin loop corrections, 
and denoted by m^. In this paper, we take 



m„ = md = m, = 0, = 1.61^^2 GeV, = 4.791^;^^ GeV. 



:i+o.i9 



(37) 



As for the B-meson lifetimes, we take [H] tb^ = 1.638 ps and tb^ = 1.530 ps, respectively. 



C3. The decay constants and form factors 

In this paper, we take the heavy-to-light transition form factors [51] 

F(f-^^(0) = 0.258 ± 0.031, F^^^'{0) = 0.331 ± 0.041, V^^^'{0) = 0.411 ± 0.033, 
A^^^* (0) = 0.374 ± 0.034, Af^^* (0) = 0.292 ± 0.028, ¥^^''{0) = 0.323 ± 0.030, 
^^^''(0) = 0.303 ± 0.029, Af^^(O) = 0.242 ± 0.023. (38) 

and the decay constants 

/b = (216 ± 22) MeV [50] , = (130.4 ± 0.2) MeV [19] , /i^ = (155.5 ± 0.8) MeV [19] , 
fx* = (217 ± 5) MeV [51], fp = (209 ± 2) MeV [51]. (39) 



C4. The LCDAs of mesons and Ught-cone projector operators. 

The light-cone projector operators of light mesons in momentum space read [521 [3] 



ifp 
4 



fc2 ■ ki 



a(3 



{M() 



ifv 



mvfv 



fv h ■ h 



(40) 



(41) 



af3 



where fpy are the decay constants, and fip = mbr^/2, with the chirally-enhanced factor 



defined as 



2ml 



mf,(/i)2mg(/i) 



2m|. 



(42) 
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where the quark masses are all running masses defined in the MS scheme. For the LCD As of 
mesons, we use their asymptotic forms [531 El! 

^py{x) = 6x{l-x), Mx) = l, Mx) = H'2x-l). (43) 

As for the B-meson wave function, we take the form [55] 

^BiO = Nb^I - Oexp[ - (xr^)'(e - ^bY] , (44) 

L \Mb — TTli, / J 

where = 1 — 'mb/Ms, and A''^ is the normalization constant to insure that d^^siO = 1- 
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